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Abstract 

We consider boundary Harnack inequalities for regional fractional Laplacian which are 
"q«^' ■ generators of censored stable-like processes on G taking 

(N ■ _ _ 

I ' as the jumping measure. When G is a C^'^ ^ open set, 1 < a < /? < 2 and k S C^(G x G) 

^ . bounded between two positive numbers, we prove a boundary Harnack inequality giving 

' ' ' dist{x,dG)'^~^ order decay for harmonic functions near the boundary. For a C^'^~^ open 

. set D C D C G, 0<a< (IV a) </3<2, we prove a boundary Harnack inequality giving 

' ' dist{x, dD)°'^^ order decay for harmonic functions near the boundary. These inequalities 
are generalizations of the known results for the homogeneous case on C^'^ open sets. We 

' also prove the boundary Harnack inequality for regional fractional Laplacian on Lipschitz 

p I . domain. 

■ Key words fractional Laplacian, symmetric a-stable processes, censored stable processes, 

I boundary Harnack inequality 

g '■ MR(2000) Subject Classification: Primary 31A20, Secondary 60G52, 60J75 

[ 1 Introduction 

>: 

. Let G be an open set in and k a positive symmetric function on G x G. For < a < 2, the 

I regional fractional (fractional-like) Laplacian is defined by 

in ■ AA' u(x) = lim^(n,-a / ^ '7^ ^ ^ dy, x e G, 1.1) 

i> : 

O ■ provided the limit exists, see [25]. Here A{n,-a) = |a|2""^r((n + a)/2)7r-"/Vr(l - a/2) 

' coming from = — (— A)"/^ when k = 1 and we refer to Guan and Ma [24] for k = 1 

^ ■ in (jl.ip . Under some regularity conditions, it is known that the a/2 power of a second order 

elliptic operator with Neumann boundary condition is an example of (II. ip . Since the integral 
kernel in (jl.ip may not be homogeneous in space, these operators to fractional Laplacian are 

similar to the second order elliptic operators to Laplacian. For 1 < a < 2, among others, an 

-,i 

explicit boundary Harnack inequality (BHI) for A-^' was proved in Bogdan, Burdzy and Chen 
[12] on C^'^ open sets, where it is called the BHI of the censored stable processes. The main 
aim of this paper is to consider the same type inequality for the nonhomogeneous case and the 
corresponding BHI on Lipschitz domain. 

Boundary Harnack inequalities are important tools in studying the boundary value problems 
in partial differential equations and potential theory of Markov processes. Analytically, such 
inequalities describe an uniform asymptotic behavior for solutions of the Dirichlet problems 
near the boundary. In Chen and Kim [T3], the BHI in |12j was used in the proof of the Green 
function estimates of censored stable processes. See also Bogdan [9] for the Brownian motion 
case. We refer to Bass |2j, Chen, Kim and Song [15] [16] for more applications of the BHI. 
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Boundary Harnack inequalities were first proved for Laplacian in Dahlberg [19j and Ancona 
[1] on Lipschitz domains. It was later extended to the second order elliptic operators in diver- 
gence form in Caffarelli, Fabes, Mortola and Salsa [13], and in nondivergence form in Fabes, 
Garofalo, Marfn-Malave and Salsa [20j . A probability method for studying such inequalities 
began in Bass and Burdzy f^. This method was applied to prove the BHI for Laplacian on 
Holder domains for elliptic operators in divergence form in Banuelos, Bass and Burdzy 

The study of the BHI for the fractional Laplacian began in Bogdan [8], Bogdan and By- 
czkowski |10) on Lipschitz open sets. Significant progresses have been made on open sets in 
Song and Wu [28] and the recent paper Bogdan, Kulczycki and Kwasnicki [11]. An explicit 
BHI for the fractional Laplacian was first given in Chen and Song [I8j on C^'^ open sets. Due 
to the jumps of stable processes or equivalently the nonlocal property of their generators, the 
corresponding harmonic functions show different feathers from the Laplacian case. Compared 
with the fractional Laplacian case in |18] (0 < a < 2), the BHI in [12] for the regional frac- 
tional Laplacian (1 < a < 2) gives a different decay for harmonic functions near the boundary, 
i.e., the former is of order p(x)°'^'^ and the later is of order In [12], the Markov pro- 

cesses associated with the regional fractional Laplacian under the Dirichlet boundary condition 
were first introduced and called the censored stable processes. We refer to [12] for some other 
characterizations of this process. 

A standard box method to prove the BHI includes comparison of harmonic measures around 
boxes, the Harnack inequality and the Carleson estimate. We refer to Bass and Burdzy [3] for 
this method in the diffusion case. The proof of the BHI in [12j studied these steps mainly by 
explicit harmonic functions given in the same paper and a relation between the censored stable 
processes and the symmetric a-stable processes. Some strong techniques are involved in this 
original proof. Due to the importance of this inequality, it is helpful to simplify the proof in [12] 
and to study this result in more general situations. In particular, our arguments can be used to 
study the Lipschitz case which is an open problem in this direction. We remark that the (super 
sub) harmonic functions given in [12j plays a fundamental role in this paper. 

To further introduce the results and the methods of this paper, we prepare some definitions 
below. For x = (xi, • • • , x„) G M", we write x = (x, x„). Let < 7 < 1 and F : M"^^ M. We 
say that F is a C^''*' function if it is differentiable and 

|VF(y) - VF(£)| , , 

||F||i,^:= sup J 1~ ^' <oo, (1.2) 

y^x,\y-x\<2 W — X\' 

where V = {d/dxi)^~^. The constant 2 in (|1.2p is only for the convenience of the later use. 
Let G be an open set in M". We say that G is a special C^'"^ domain if for some C^'"^ function 
F : R"-"^ ^ M, G can be represented as {x = {x,Xn) G W^,Xn > F(x)}. In this case G is also 
denoted by Gr- We say that G is C^'^~^ if there exist vq > and A > such that for each 
z G dG, we can find a C^'^~^ function F^ : — > M with ||F||i^^ < A and an orthonormal 
coordinate system GSz such that 

GnB{z,ro) = {y = {yi,--- ,yn) ■■ y„ > F^(yi, • • • y„_i) } n ^(z, tq) . (1.3) 

By rotation and translation, we can always assume that WzCz) =Tz{'z) = in CSz- The pair 
(ro, A) is called the characteristics of G. The characteristics of a Lipschitz open set is defined in 
a similar way. For each 6 > 0, set 

G's = {yGG: p{y) <5], Gs = {y^G: p{y) > 6}, (1.4) 

where p{y) = dist{y,dG). 

Let 1 < a < 2 and 'ipi,'4'2 be positive functions in G^{G x G). Let k be a symmetric function 
on G X G taking values between two positive numbers Ci and C2. Assume also that for some 
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constant C" > and 6 G (0, tq) 



K{x,y) - iJi{x,y) - 4^2{x,y) \l_l\n+c. <C'\x-y\, x,y£G'g, 
[ \n{x,y) - k{x,x)\ <C'\x -y\, x,y e Gs/2, 

where y is the reflection point of y with respect to dG (see section 4). Write 

M:=G'+ sup {\VyMx,y)\ + \^yMx,y)\) 

x,y£G,\x—y\<l 

and denote by (Xt) the reflected stable-hke process. The following theorem is an extension of 
the BHI in |12] for k = 1 on G^'^ domain. 



Theorem 1.1. Assume that a,K satisfy all the conditions above and G is a G^'^ open set in 
M" with characteristics < 1 and A. Let Q G dG and r E (0,ro). Assume that u > is a 
function on G which is not identical to zero, harmonic on Gr\B{Q,r) and vanishes continuously 
on dG n B{Q,r). Then there is a constant G = G{n,a,A,6,Gi,G2,M) such that 

^<C4^, x,yeGnB{Q,r/2). (1.6) 

Hv) p{yr 

Moreover, if ip2 = Q in this boundary Harnack inequality holds for G^'^~^ open sets with 

l<a< p <2. 

Here the notation G = G{n,a,A,5i,Ci,G2,M) means that the constant G is positive and 
depends only on parameters in the bracket. This convention will be used throughout the paper. 
When ■02 = in (|1.6|) . the last conclusion in Theorem 1.1 was conjectured in [T^]. We remark 
that P = a is the critical value in our proof and Theorem 1.1 may not hold for this value. 

In [12], when G is a special G^'^ domain and k = 1, to get sharp estimates of harmonic 
measures, (super) subharmonic functions are constructed by explicit harmonic functions on M!j: 
and non-explicit perturbations deflned through the symmetric a-stable process on M". Here we 
construct explicit (super) subharmonic functions by perturbation directly. This construction 
may also be used to prove the known explicit BHI for Laplacian, i.e., a = /3 = 2 in (jl.6p . 

For the Harnack inequality, we may adopt the method in Bass and Levin [6]. Here we give 
another proof which might be more straightforward for these nonlocal operators. This proof 
is similar to the proof of the Carleson estimate given in Lemma 14.21 which is an application of 
the box method for jump processes taking (j4.16p as the key observation . We remark that the 
method in [6] can be applied to prove the Harnack inequality for jump diffusions, see Song and 
Vondracek [30]. Therefore by applying the method in this paper, we may prove the BHI for 
operators like A + A°/^ on G^'^ open sets, where the decay is of order p{x) near the boundary. 

The BHI for the fractional Laplacian on C^'^ open sets was proved by Poisson kernel estimates 
in [18]. This and many other estimates of the symmetric stable processes given before depend 
on their explicit Poisson kernel and Green function for a unit ball which are not available for 
the nonhomogeneous case. In Lemma l6. 11 we present explicit (super, sub) harmonic functions 
of the fractional Laplacian on half spaces which allow us to study the nonhomogeneous case. 
See Theorem 16.41 We notice that the harmonic function in Lemma 16.11 has been obtained in 
Banuelos and Bogdan [7]. As applications, for the fractional-like Laplacian under condition 
(j6.7p , we may get the sharp estimates of their Green function and Poisson kernel as in [18] and 
hence we may get the BHI in [llj under the same conditions. 

Another main result of this paper is the following boundary Harnack inequality on Lipschitz 
domain. The strategy of the proof is essentially the same as the proof of Theorem 1.1. 
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Theorem 1.2. Let G be a Lipschitz open set in M" with characteristics rg < 1 and A. Assume 
that 1 < a < 2 and k be a C^{G x G) function bounded between two positive numbers. Let 
Q € dG and r E (0,ro). Then there is a constant C such that 

^<C^, x,y^Gr^B{Q,r/2), (1.7) 
u{y) v{y) 

where u,v >0 are functions on G which is not identical to zero, harmonic on G H B{Q,r) and 
vanishes continuously on dGriB{Q,r). 

To prove Theorem 1.2, the heat kernel estimate of the reflected stable processes in Chen 
and Kumagai [T7] is used to give some hitting probability estimate. The censored stable pro- 
cesses can be extended to the reflected processes on G which is formulated in [12]. For general 
AC, these two processes are called the censored stable-like process and the reflected stable-like 
process respectively (see Remark 2.4 [12j). They are symmetric Markov processes on G and G, 
respectively. The Dirichlet form of the reflected stable-like process is 

^ ^ 2 ^ J JgxG 

^« = {uGL2(G): £'^{u,u) <oo}, (1.8) 

where k{x, y) is bounded between two positive numbers and G is Lipschitz. In [1^, this reflected 
process was refined to be a Feller process (Xt) on G under a more general condition. 

The structure of this paper is the following. In Section 2 we study (super, sub) harmonic 
functions. In Section 3 we prove the Harnack inequality. In Sections 4 and 6 we prove boundary 
Harnack inequalities for the regional fractional-like Laplacian and the fractional-like Laplacian, 
respectively. The Lipschitz case is studied in section 5. For a, 6 G M, a V 6 := max{a, 6}. We 
use m(-) to denote the area measure of (n — l)-dimensional subset. For any set U, denote 
Tjj = inf{t > : Xt ^ U}. For function u on M", we always take it as a function on G by 
restriction when considering A'q'^u. The dimension n is assumed bigger than two throughout 
the paper. 

2 (Super) subharmonic functions for regional fractional Laplacian, k = 1 

Let u he a Borel function on G and let U be an open subset of G. We say that u is a (super, 

sub) harmonic function on U with respect to A^' if A^,' u{x){<, >) = for x ^ U. We say 
that n is a (super, sub) harmonic function on U with respect to the reflected stable-like process 
m if 

u{x){>,<) = EMXrs), xeB (2.1) 

for any bounded open set B with B C U, where tb = mi{t > : Xt ^ B}. Under the conditions 
of K and G in this paper, the harmonic function for the reflected stable-like process is continuous 
on U (see Corollarv 13.61 below). This implies that it is harmonic for Aq in the weak sense 
(cf. Theorem 6.6 [25j for k = 1). When u is and k is C^, Theorem 4.8 in [23j shows that 
these two definitions are equivalent. In what follows we use (j2.ip for the definition of harmonic 

functions. We write A^' by A^ when x £ G. 

In [12], to establish the BHI for k = 1, the following estimates are given for the regional 
fractional Laplacian acting on function u = p"~^: 



A''fuix)<Apixf-^ x€G[/^, iff3<2, 
\A%^\{x)\<A\lnp{x)l xeG[^^, if 13 = 2, 



(2.2) 
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where G is a special C^'^ ^ domain and A is a positive constant (1 < a < /? < 2). When /3 < 2, 
we can not get a bound for A^'^^u because it may take — oo. This is related to the fact that p 
may not be when dG is C^'^~^ (/? < 2). To improve estimates ()2.2p . we replace the distance 
function by a "height function" which is equal to x„ — Tz{x) in a neighborhood of z E 9G. 

Since there is no difference for boundary conditions discussed here when n = 1, we always 
assume that n > 2. Let Wp{y) = Un for y G R" and p £ M. Our starting point is the following 
explicit harmonic functions and ()2.16p below given in (5.4) |12| 



A ' Wa-l[x) 



0, X e 



a G (1,2). 



(2.3) 



In the following lemma, when necessary, a function defined on a domain is also considered as a 
function on M" by taking zero outside. 

Lemma 2.1. Let 1 < a < f3 < 2 and let T : M""! be a C^'^-i function with r(0) = and 
Vr(0) = 0. Define function h^-iix) = {x^ — r(x))"~"^/||j|<2} /o?" x £ D := D^- Then there 
exists constant Ai = ^i(n,a,/5, |[r||i^^_i) such that 



|A^/'/i,_i(x)I < 



Aip{xf-\ 
^i(Ilnp(x)[ + l), 



X G D[, \x\ < 1, if (3 <2, 
X G D[, \x\ < 1, if 13 = 2. 



(2.4) 



Proof Denote ha-i by /i. We only prove the lemma for a < (3 < 2 because the proof for 
/? = 2 is similar. Let x G D[ with \x\ < 1 and choose a point xq G dD satisfying x = xq. 
Denote by li(xo) the inward unit normal vector at xq for dD and set <&(y) = {y — xo,T?(xo)) 
for y G M". It is clear that 11 = {?/ : $(y) = 0} is the plane which is tangent to dD at point xq. 
Let r* : X G M"~^ ^ M be the function of plane 11, i.e., 



and set 



((x, r*(x)) - xq, r?(xo)) = 0, 



U = {y=iy,yn):y€D, \y - x\ < 1, y„ < 2 + 2^||r|li,^_i }. 



Write h{y) = |y„ — T*{y)\ for y G M". Applying the assumption that dG is G^'^ ^ and Vr(x) 
Vr*(x) = 0, we have by the mean value theorem 



\hiy) - h—{y)\ < |r(y) - T*iy)\ < ||r||i,^_i[y 



1/3 



2/G [/. 



(2.5) 



Let p ^{y) = dist{y ,U) for y G M" and Dr* 
h = y/1 + [Vr(SJ5)[2pn- So we have by ([23]) 



A^/i" '(y) = (i + |vr(55) 



= {y G M" : y„ > r*(y)}. It is clear that 
'-^Al^^p^-\y) = 0, yeDr*. (2.6) 



Denote 



A = {y: T*{y) < y„ < r(y), |y - x| < 1} U {y : T{y) < y„ < r*(y), |y - x| < 1}. 

h{x) and B(x, 1) n C f/ (by the fact that x„ < 1 + 2'3||r||i_^_i), we 



Noticing that /i" ^(x 



have by (j2.6 



lim 



h{y) - h{x) 



'y£D,\y-x\>e \^ V 

-a-1 

/ / 

< lim sup / 

£j,0 JyeU,\y-x\>£ 



n+a 



dy 



r-\x) 



x-y 



\n+a 



dy 



+ lim sup 



Ky) 



r-\y) 



y&U,\y-x\>e P y 



\n+a 



dy 
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lim sup 

elO 



I 



h{y) - h{x) 

y(.D\U,\y-x\>s \^ - 2/h+" 



dy 



< 



\x - y\ 



n+a 



dy + 



\r-\y)-r-\x) 



\h{y) - h'^-'iy)] 
\x - y\"-+°' 

:=/i + l2 + h+ h- 



dy + 



B(x,iY \x-y\ 

\h{y)-h{x)\ 



n+a 



dy 



B(x,iY 1^ ~ y\ 



n+a 



dy 



(2.7) 



Noticing that Ad {y ■.\yn - (a:o)n| < 2^ ^||r||i,/3_i|y - xq\ }, we have 

\x-y\>{l + 2^^-^\\T\\lp_^)-y^h{x) > (1 + 2/3"i||r||i,^„i)-iMx), y G A, 



which imphes 



I , ^ (1 + 2/^-1 ||r||i,;3_i)-i/^(x) + |^-£| 
\x-y\> , y ^ A. 



By ()2.5p . we also have h{y) < ||r||i^^_i|y — x\^ for y & A. Therefore 



h < dr 



J\y~x\=T 



<(l|r|I?/i + i) 



\x - y\ 



h{x)P M 



dr 



\y-x\=r 



hixY^ 



(||r||"-i, + 1) _ , dr 



r./3(Q-l) 



h{x)V M J\y-x\=r 1^ Vl 



n+a 



m{dy) 



<(2^r(||r||°^_i + ||r||i,^_i) 
+(2^r(lir||Vi + lir||i,;3-i) 



~^^''^^''\^^^a-i^ r + {l + 2^-^\\ni,p_,y^h{x) ^^^ 



lh{x)'^Al 

^ tn„\n n\T^\\a , iit^ii \r)a~f3+2i'^~^'^ 



,/3(a-l),^ + (l + 2^-ir||i,;3-l)-^/i(x),_. 



+0-2 

+^-2 dr 



\a-/3+l 



a- (3 + 1 



h(x 



1/3-2 



+(27rr(||r|jVi + l|r||i,/5-i)2""^+' /_ . r 



al3~a-2 



dr. 



h{x)VAl 

As a/3 - a - 2 > /32 - 2/3 - 1 for 1< a, /3 < 2, we get 



(2.8) 



Jh(x)^ 



,al3-a-2 



dr < 



h{x)^Al 



^ {h{x) A 1)^-2. 



2/?-/3- 



//i(a;)^Al 

The following properties follows from the definitions of h and h. 

IMx) - %)| < (1 + ||r||i,^_i)|x - y|, ye M", 

|/i(x)^ -/i(y)^| < (l + 2^-i||r||i,^_i)|x-yI, y G L», |y| < 2, 

P(y)"-' < h{y) < (1 + 2^-i||r||i,^_i)"-V(y)"-\ y g \y\ < i. 

Noticing that p{x) < 1 and h{y) = for |y| > 2, by (l2T0]l - (l2T2|) 

/2 + /4 



(2.9) 



(2.10) 
(2.11) 
(2.12) 
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(i + ||r||i,^_i)"-i ^ r (i + 2/^-ir||i.^_ir-i ^ r h{x) ^ 

~ Jb(x,iy - yl"''' ^ Jb{'.,iy \x - ^ Jb{x,iy - ^ 

<(27rr ((1 + ||r||i,^_i)"-i + (1 + 2^-i||r||i,^_i)-i + (1 + 2^-i||r||i,^.i)"-i). (2.13) 

To estimate Is we define a transform '^{y) = (J, Zn) by 

z = y, Zn = yn-r*{y), yGM". 
We see that |^| = 1, where ^ is the Jacobian determinant of ^. We can also check that 

\yi - y2\ > (1 + ||r||i,^_i)-i|^(yi) - ^-(^2)1, for yi,y2 G M", 
and ^{U) C {2/ : |y - x| < 1, |y„| < 2^+i(||r||i^^_i + 1)}. Hence by i^E), the inequahty 

^ ^a-2|^_^|^ 6 > 0,0 > 0,1< a < 2, (2.14) 
and applying the transform ^, we have 

. < /• (i + iiriK,^_ir+°+i|%)r-^|y-g|/^ 

(i + ||r||i,^„i)"+°+^|z„r2|?-i|/^ 

*{[/) |(£,7^(x)) - 2;["+" 

.2/3+i{||r|K„,_i+l) . _^^n+a+l|^|a-2|j_~|/3 

< / dr / _ m[dz) 

J-2/3+i{||r||i,^_i+i) J zr,=T,\z-x\<\r-h{x)\ |(a:,/i(x)) - 

.2^+H\\ni,,-i+l) f (i + ||r|| _^)n+a+l|^|a-2 

+ / or / _^ m(d2;) 

J_2/3+i(||r||i,^_i+l) J z„=r,\r-h{x)\<\I-x\<l \z - x\'^+°' P 

^ (27r)" /•2'^+MI|r||i,,-i+i) (i + |[r||i,^_ir+"+i[r|"-^ 



(27r)" /•2'3+^{||r||i,,-i+l) (-^^|[p|[^^^_^)n+a+l|^|a-2 



J IZlli / ^ ^' LJ J. 

a-Z^ + Uo |/i(x) -r|"-/5+i 

, 2(2.r(i + ||r||,,_,r"+^ / %)^-^ ^Tlfx)/^-^ (2 15) 

Combining (fO) - (12:9]) . (12131) and (l2l^ . we get ([23)). □ 

Remark 2.1. Estimates 1^2.4^ may not hold if we take (3 = a in Lemma \2.1\ For n = 2, 

r(xi) = Ixil'^ and x* = (0,t) with t > 0, we can check that dy = —oo and 

-^1,-^2,-^4 (i-re all finite. This gives A'^^hix*) = -oo. When a < (3 < 2, we can also prove that 

Al)^p'^~^{x*) may take — oo. We still consider the above example. Let Xq he the point on dDr 
such that \xq — x* \ = p{x) and (xq)i > 0. Let 

U = {{yi,y2) ■■ y2 > \yif or yi < 0} n {(yi,y2) : y2 > 0} 
and denote the distance function to dU by £,{x). Since is smooth in a neighborhood of x* , we 
know that A^^"~^{x*) is finite. On the other hand, 

r pjyr-'-qyr-' 

JDr 1 3 

a. 

Hence we have A^^p"~-'^(x*) = — oo. 



-^2+a dy = -OO. 
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Recall Wp{y) = yl for y G W^. By (5.4) in [E] 

AS{,'^x;p(x)=^(n,-a)^^S(^,^)7(a,pK-", x € M!^, pG(-l,a), (2.16) 

where Wn-i is the (n — 2) -dimensional Lebesgue measure of the unit sphere in M""^, B is the 

Beta function and 7(0, p) = (^x_iY^a dt. In what follows we denote the constant on 

the right hand side of (|2.16p by C{n,a,p). 

Lemma 2.2. Let a, T and D be described in Lemma \2.1\ and let p be a number such that 
a > p > a — 1. Define function hp{x) = {xn — T{x)Ylj^\^^2} on D = Dy- Then there exists 
constant A2 = A2{n,a, f3,p, ||r||i^^_i) such that 

A^/'/ip(x) > ^2p(x)P-", X G D'y^^, \x\ < 1. (2.17) 

Proof We use the definitions and the notations in the proof of Lemma 12. 1[ Following the 
arguments in ()2.7p . for x £ D[ with \x\ < 1 we have 

hm / h iy)-h {x) 

slO Jy^D,\y-x\>e \x - 2/1"+'" 

^/ N r, rv f \h^ iv) —'hf' {x)\ , f \h^ (y) — h^ (x)\ , 

>C(n,a,p)xP-°' - / — —^dy- / — , ,^ dy 

J A \x — y\ ^ Jb{x,iy f — 2/1 



hp{y)-h{y)\ f \hp{y)-hp{x) 



,n+. - / ' 7^ i^^a ^2/ 

[7 |x-y|"+" Jb(x,i)- 

=C(n, a,p)xP-" -h-h-h- h. (2.18) 
By similar calculations as in Lemma 12.1^ we can find constant ki such that 
h < A:i(p(x)'^+P-"-i V 1 + I Inp(x)l), 

h < fci(p(x)'^+P-"-i V 1 + I Inp(x)l), l2+h< h. (2.19) 
Noticing that p - a < 0, /3 > 1 and C(n, a,p) > 0, we obtain ([2171) by (l2J8]) and (l2J9]) . □ 

Lemma 2.3. Let a, T, D, h^-i and hp be objects described in Lemmas \2.1\ and \2.2l Let f be a 
bounded function in C^{D). Then there exists constant A3 = ^3(71., a,p, sup \ f{y)\, sup |V/(y)|) 

j/eM" \y\<2 

such that for X e Dn B{0, 1) 

\{fiy)-fix)){hp{y)-hpix))\ f ^3(|logp(x)| + l), p = a-l, 

^ - \ ^3, a-l<p<a. ^' ' 

Proof We only prove the lemma for p = a — 1 because the others can be proved similarly. 
Denote h^-i by h and let x G L> n B{0, 1). By (|2TT|) . (12^21) and (imH . 



\if{y)-f{x)){h{y)-hix))\ 



1^.. M f \h(y)^ - h{x)^f-^ , 
< sup |V/(y)| / |„Vli ' dy 

\y\<2 J Dr){p{x)<\y-x\<l} W - X\ ^ 

1 1 a-2 

+ sup |V/(y)l / LML^ ^ ^ J dy 

\y\<2 JDnBix,p{x)) 
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+2 sup \f{y)\ / }^^^ dy 

JDnB{x,l)<'n{y:\y\<2} W — X\ 



+2 sup \fiy)\ I dy 

yGR" J{y:\y\>2} \y " X\^^°' 

< sup |V/(y)| / ■ dy 

\y\<2 J Dn{p{x)<\y-x\<l} lU ~ X\ 

+ sup |V/(y)l / ■ -—- dy 

\y\<2 JDnB{x,p{x)) 

+2 sup \f{y)\ / ^ —: dy + 2 sup \f{y)\ {2tt) 

S/GK" JDr\B{x,lY \y — ^ j/GK" « 

<(sup |V/(y)| + 2 sup + 2^||r,||i,;3-i)"-'(27r)"(- Inp(x) + + 1 + -), 

\y\<2 2/GK" ^ - a a 

which completes the proof. □ 
For function k on G x G and Q G in the proposition below we denote 

Co= sup \VyK{x,y)\. (2.21) 
a;,j/G-B((9,i)nG 

Proposition 2.4. Let 1 < a < /3 < 2 and G a C^'l^~^ open set in with characteristics tq = 1 
and A. Let k be a function on G x G taking values between two positive numbers Ci and G2. 
Then for a — l<p<a and Q G dG, there exist function Up and positive constants A4 = 744(A), 
= (n, a, j3,p, A, Go, Gi, G2) such that 



^4 ^lGnB{Q,2/?,)P{xY < Up{x) < A4lGnB{Q,2/3)P{xT, x G G, (2.22) 



and 



A^' Up{x)>A5p{x)P-'^, X eGnB{Q,l/A5), a - 1< p < a, (2.23) 

(r)\<i ^5P(x)^-^ xGGni?(Q,l/2), i//5<2, 
I^G ^1 ^5|lnp(x)|, xGGn5(Q,l/2), ifP = 2. ^^'^^^ 

Proof Without loss of generality, we assume that Q = and take the coordinate system CSq (see 
([L3D). Define functions lip (x) = ixn-TQix))PlGri{B{Q,2/3)} onG andvp{x) = (x„-rQ(x))P/[j|<2 
on for a — 1 < p < a. It is easy to see that (|2.22|) holds. When k = 1, noticing that for 
X G G n -B(0, 1/2) the integral in 1^ for on G n B{Q, 2/3)^ and Vp on Dr^ n B{Q, 2/3)^ can 
be bounded by constants depending on n and q, we can prove this proposition by Lemma 2.1 
and Lemma |2.2[ For general cases, the conclusion can be proved by the case k = 1, Lemma 2.3 
and the following identity: 

a|'^/i(x) 

f (t^(x,v) — K(x,x))(h(y) — h(x)) , , ^.nl2,, s 
=Ain,-a)lim ^ ^ , ^-^dy + k{x,x)A'^^\{x). (2.25) 

^i-^ JyGG,\y-x\>e F ~ Vl 

□ 

— K 

3 Harnack inequalities of A^' 

The following example can be found in 
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Example 3.1. LetG = Wlandy = {y,-yn)fory = {y,yn). For k{x, y) = 1+ ; A,^,' is 

the formal generator of the subordinate reflected Brownian motion on . When G = (0, 1) and 

K{x,y) = Yl'k^=-oo 1-^ ~ yl^'''"/!^^ i y + 2k\^'^°', A^' is the formal generator of the subordinate 
reflected Brownian motion on [0, 1] . 

Remark 3.1. Define function Wp{y) = yn for y € M". When K{x,y) = ji^_|jn+a > we have 
fseel28^) 

A|ru;p(x)=^(n,-a)i^^(^,^)7(a,|.K-", x G M^, j>G(-l,a), (3.1) 

where j(a,p) = (t^ — 1)(1 — t"~P~^)/(l + t)^~^°' dt. This gives the same (super, sub) harmonic 
functions as the homogeneous case in i2.16\) . which will be used later. 

Notice that the derivatives of k, in the examples above are not bounded. To give results 
including these examples we introduce the following condition. Let < Ci < C2, C3 > and 
7 < 0. We say that k or the reflected stable-like process (Xt) satisfies condition [Ci, G2, (^3,7] if 

Ci < K{x,y) < C2, x,yeG; y) - k(x, x)| < C3(p(x)'^ V l)|x - y|, x,y£G. (3.2) 

We can check that functions k in the Example 3.1 above satisfy condition [Ci, C2, C3, — 1] for 
some constants Ci, C2, C3 > 0. 

Next we prepare a stochastic calculus formula for (Xt). For a measurable function / on G, 
denote / G £i{G) if 



sup 



/ 7-. ] — ; — dy < 00. (3.3) 

(l + Ix-2/|)"+" ^ ^ 



For any subset U C and < 7 < 1, we say that u is uniformly 7-Holder continuous on U if 

\u(y) — u(z)\ , 
sup ^-f^ < 00. (3.4) 

We shall denote {u E G^'^{U)) u G C^{U) if (all the first derivatives of u) u is uniformly 7- 
Holder continuous on U. For any 6 > and A C M", define ta = inf{t > : Xt A'^} and 
= {y ■ \y — x\ < ^5 for some x € A}. For any relatively open subset A of G, we denote 
by (pf ) and G^ the probability transition function and the Green function of (Xt) killed upon 
leaving A, respectively. In [23], a semi-martingale decomposition of f{Xt) is given for / € Gl{G) 
(see |24) for the homogeneous case). To consider more general functions, we prove the following 
results. 

Proposition 3.1. Let G be a Lipschitz open set in M". For 1 < a < /? < 2, let k be a symmetric 
function on G x G satisfying condition [Ci, C2, C3, 7] with > 7 > a — 3. For < a < (3 < 1, 
let K be a measurable symmetric function on G x G bounded between Gi and C2. Then for f 
belonging to 

n Ci{G), l<a<p<2- C'^i'G) n Ci{G), < a < /3 < 1, (3.5) 

we have 

f{Xt) = /(xo) + Mt+ f Al;'''f{Xs) ds, a.s. xq S G, (3.6) 

Jo 
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where {Mt)t>o is a martingale. If A is a relatively open set in G and, for some 5 > 0, f satisfies 
US. 5\) with 'G replaced byGCiA^, then 



E.^ifiXt^r^)) = fixo) + E,,i A^' f{Xs)ds), xoeA,t>0. (3.7) 

Jo 

Moreover, if Pxq{Xt-^ £ dA) = and f is a positive function such that / = on A, then (3.1) 
still holds. 



Proof Assume that / satisfies (I3.5p . For 1 < a < /3 < 2, by (j3.4p . the derivatives of / at point 
X G G can be bounded by ki{l + for some constant ki > 0. By this estimate, (j3.3p and 
straightforward calculations for the integral in (jl.ip on sets B{x,p{x)), {GnB{x, l))\B{x, p{x)) 
and G \ B{x, 1) respectively, we can prove that for some constant k2, 

|Aj'V(a;)| < k2{l + |x|)p(2;)(2+T-")^(i-"), x e G[, l<a<P<2. (3.8) 
By ()3.3p and ()3.4p we can find constant ^3 such that 

|Aj'''/(a;)| < A:3, x e Gi, 1 < a < /3 < 2. (3.9) 

By calculating the integral in (jl.ip on sets G D B{x, 1) and G \ B{x, 1) respectively, we can also 
check 

|Aj'''/(x)| </c4, x£G, 0<a<(3<l, (3.10) 
for some constant k^. Noticing that (2 + 7 — a)A(l — a) > —1, with the help of the heat kernel 
estimates in |T7] and (|3.8p - (j3.1Up we can prove that -E'z(/o fi-^t)\d^) is a bounded function 

on G (c.f. Lemma 4.6 p^). This implies that Ex{Jq \Aq f{Xt)\dt) is a bounded function 
on G for any t > 0. Thus we can prove (j3.6p by Theorem 5.25 [21] at time t A T^(o,n) (c-f- 
Theorem 4.1 |24j) and letting n — > 00. Formula (13. 7p is a consequence of ()3.6p by approximation 
procedure. □ 

For a relatively open set A in G, we say that A has outer cone property in G if, for some 
r] > and each Q ^ dA, there is a cone in G\A isometric to {x E M" : \{xi, ■ ■ ■ ,Xn-i)\ < vl^nl} 
and taking Q as the vertex. 

Proposition 3.2. Let a,G and k be the same as in Proposition \3.1\ Let A Q G be a relatively 
open set with outer cone property in G and define r = mi{t > : E A^}. Then the 
distribution of X^ is absolutely continuous on G\A when (Xt) starting from A. For any t > 0, 
we have 

Rx{XrI{r<t} e dy}/dy 
=Ain,-a)J^^J^^^^^0^^dzys, {x,y)eAx(G\A). (3.11) 

Furthermore, 

P,{Xr E dy}/dy 

=A{n, -a) [ ^i^Ml^l^gildz, (x, y)eAx(G\A). (3.12) 

Ja \^ y\ 

Proof To show that Px{Xt- E OA} = for x E A, by the method in Lemma 6 [8|, we only need 
to prove that there exists a constant c such that Px{tb{x,p{x)) E G \ A} > c for any x E A (the 
boundedness assumption in [8] is not necessary because A can be approximated by A n B{0,n) 
by letting n | 00). We omit the proof of this estimate because it is similar to (j3.25p below. Thus 
we can prove (I3.1ip by Proposition 13.11 Formula (I3.12p is a consequence of (13. lip . □ 
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Lemma 3.3. Let < a < 2 and let G be a Lipschitz open set in M". Assume that k is a 
symmetric function on G x G satisfying condition [Ci, G2, C3, —1]. Let A > 1 and define process 
{{Zt)t>o,Qxo) = {{^Xx-at)t>o, Pxo/x) for xq G AG. Then (Zt) is a reflected stable-like process 
on XG satisfying condition [Ci, C2, C3, — 1]. 

Proof The conclusion can be proved by checking that the jumping measure of (Zt) is 

^P^dxdy, x,y^XG. □ 
\x — y\^+°' 

Lemma 3.4. Let < a < 2 and let G be a Lipschitz open set in M" with characteristics = 1 
and A. Assume that n is a symmetric function on G x G satisfying condition [Ci,C2,C3, — 1]. 
Then for < e < 1, there exists constants Aq = AQ{n,a,C2,C3,£) and A'^ = ^Q(n,a,Ci) such 
that for any < r < ro/2 

A%f°' < inf EyT^/^j.\ < sup EyT^ix r) < ^6'^°) X G G with p{x) > 2r. (3.13) 

yeB(x,(l-e)r) ' 

Moveover, the last inequality in hS.lS^ holds for all x ^ G provided r < r^/A, where A'q depends 
further on A. 

Proof By the scahng property in Lemma 13.31 and the Lipschitz condition of G, we can assume 
that r = 1. Choose /i G C'^iG) such that < /i < 1 and 

/i(y) = 0, y(^B{x,l-e)- = 1, y ^ B{x,l - e/2r . 

By direct calculation, we can find a constant ki = /ci(n, a, G2, G3, e) such that |A^' /i(y)| < ki 
for y G B{x, 1). Thus we can prove the first inequality in (I3.13P by applying formula (13. 7p to 
£'y(/i(XT-g(^ )). Similarly, with the help of Proposition 13. 2|, the last inequality in (|3.13p can 

— K 

be proved by considering function /2 = Ig\b{x1)^ where we can check that A^' /2(y) > ^2, 
y G G n B(x, 1), for some constant /c2 = k2{n, a, Gi). For the last conclusion, G PI B{x, 1) may 
not have the outer cone property in G, where we need to replace B{x, 1) by a bigger set in 
B{x, 2) satisfying this property. □ 

The next theorem extends the Harnack inequality for the censored stable process in jl2) . 

Theorem 3.5. Let < a < 2 and let G be a Lipschitz open set in with characteristics tq 
and A. Assume that k is a symmetric function on G x G satisfying condition [Gi,G2,G3, — 1]. 
Let < r < 1, /c G {1, 2, . . .} and xi,X2 G G such that B{xi,r) U B{x2,r) C G and \xi — 
X2I < 2^r. If u > is harmonic for (Xt) on B{xi,r) U B{x2,r), then there exists constant 
Aj = Ai{n,a,Ci,C2,C3) such that 

^_l2-fc(n+a)^(^^^ < u{xi) < ^72^("+°)u(2;2). (3.14) 

Proof For simplicity we assume k = 1. Let y G with p{y) > r. First we prove that there 
exists a constant ki = ki{n,a) such that 

u{yi) < kiu{y2), yi,y2 G B{y,r/2), (3.15) 

provided u > is harmonic for (Xt) on B{y,r). To show this we only need to prove that 

u{yi) <kiu{y2), yi,y2 e B{y,r/2) and \yi - y2\ > r/3. (3.16) 

Approximating by functions := Ex{{u A /c)(Xt-^^^ ^j)), we can assume that u is bounded. By 
scaling we can also assume r = 1. 
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Let 2/1,2/2 S B{y, 1/2) such that \yi — 2/2! > 1/3. Suppose that u{yi) > Mu{y2) for some big 
number M and we can construct a sequence of points {xk) € -6(2/1, 1/6) such that xq = xi = 2/1 
and 

\u(xk)\>{l + 5f-^Mu{y2), |xfe-xfc_i| < 12-i(fc-l)-2, A; > 1, (3.17) 

for some (5 > 0, then the contradiction between the boundedness of u and hm^^oo u{xk) = 00 
leads to (|3.16p (here we assume that ^(2/2) > because u{y2) = imphes that n = on 
GnB(2/,l), c.f. (IXT9D bebw). 

Suppose that (j3.17p holds for k = \ and some 5, M which will be fixed later. Setting 
Bk = B{xk,2A-^k-'^) and Tk = for > 1 , we have by Proposition 13.21 and Lemma 13.41 

=A{n,-a) I I 9^^hAdzdy 



>2-("+°)^(n,-a) I G'^^ixk.z) [ 
JSk Jy 

=2-("+-)^(n,-a)(i?,,rfc) / 

Jy 



yeG\{2B,) -2/1"+" 
dy 



dzdy 



yeG\{2B,) \xk - 2/1"+° 



>ki (3.18) 

for some constant ki = ki{n,a). Similarly, by setting Bq = ^(2/2, 1/6) and tq = tbq, we also 
have 

j u(v) 

u{y2)=Ey^{u{Xra)Ix <.G\Bo)>k2{n,a) u^dy. (3.19) 

By ()2.ip and an estimate of -Pyal'^o ^ B{y, 1/2) \ {2Bq)) similar to (j3.18p . we can find 2/3 G 
B(y, 1/2) \ {2Bo) such that ^(2/3) < ^3?-t(2/2) for some constant ^3 = k3{n,a). Similar to (j3.19p . 
we have 

uiys) > h{n, a) [ ""^^l+j v- (3-20) 

JyeBo m - yr+ 

Noticing that \y — Xk\ > j^ilv — 2/2I V I2/ — 2/3I) for 2/ € G \ {2Bk), we have by Proposition 
121 and Lemma 13.41 



Ex,{u{Xr,)Ix^^eG\i2B,)) 

uiy) 



<r+^A{n,-a){E,,Tk) [ 

J V 



yeG\{2B,) \Xk - 2/1"+'^ 



dy 



dy 



yeG\{2B,) l^fc - y|"+° 

<k,{n,a)k''^( [ ^li4^dy+/ , d2/)- (3.21) 

" V2;gG\Bo - 2/1"+" Aeso 1^3 - 2/1"+" ^ ^ ^ 

By ^M-^M and ^(2/3) < A;3n(2/2), we have E^^{uiXr,)Ix,^&G\{2B,)) < kQ{n,a)k^My2)- 
Thus by dSH]), (fHTfl) and (f3l8]) we have 

(1 + 6)''-^Mu{y2) < u{xk) < (1 - A:i) sup u{y) + kek^''u{y2). (3.22) 

?;e(2Bfe)\Sfe 
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Now choose 5 = ki/2 and Kq = KQ{n, a) > 1 such that for any M > 1 

(1 + 6r~^M - hm^^ > / ~ (1 + 6r-^M, m>Ko. (3.23) 

1 — Ki/2 

If Xfc with k > Kq satisfies (I3.22P for some M > 1, then ()3.22p and (I3.23P show that there exists 
Xk+i satisfying (j3.17p for A; + 1. By (j3.22p . we can also choose M = M{n,a) > 1 big enough 
such that (|3.23p holds for 1 < A; < Kq. Therefore, we can finish the proof of (I3.15P by induction. 
Next we assume that 2'^r > \xi — X2I > r. We have for x E B{xi,r/3) 



— K f 

lB{x2,r/Z){x) =A{n,-a) I 

Jy 



yGB{x2,r/3) 1^ ~ y|"+" 



dy 



>A;7(n,a)r-"2-'=("+"). (3.24) 

By (1221), (1333]) and UjSMii we have 

^-i{^r,(.„./3) e^(^2,r/3)} 

=Exi{ A^' lB(x2,r/3){^t)dt) 

Jo 

>A;7^63~"2-'=('^+°). (3.25) 

By PJ[^ and (fX^ . 

n(xi) =i?,,(n(X,^(^^ „/3,)) > A;8(n,a)2-'=("+")n(x2), (3.26) 

which completes the proof. □ 



Corollary 3.6. Let a, G and k be the same as in Theorem \3.5\ and let u be a (Xt) harmonic 
function in an open subset D of G. Then u is continuous on D. 

Proof Let x £ D and Pd{x) = inf{|x — y\ : y € 51?}. By Theorem 13.51 we see that u is bounded 
on B(x, 2p/)(x)/3). Set r = tb[x,pd(x)/3)- ^y the strong Markov property, we have 

U{y) = EyU{Xr) = Ey[u{Xr)It>r\ + Ey[u{Xt)It<r]- (3.27) 

By the continuity of the heat kernel in [17j, we see E.[u[Xt)Ii^t<T}\ ^ ^(^(x, /Jz)(x)/3))(c.f. 
Proposition 3.6 [25]). On the other hand, 

Ey[u{Xr)I{t>r}] 

Ey[u{Xr)It>Tlx^iB(x,2pD{x)/3)\ ■ 

Therefore, by (j3.1ip . (j3.27p . Theorem 13. 51 and the dominated convergence theorem, we need only 
to check that Py{t > r} converges to zero uniformly on B{x, pd{x)/3) when t I 0. This follows 
from facts that Py{t > r} = 1 - Py{t < .} . a(i.(x, p.(^)/3)) and nmP,{. > .} = for 

y e B{x,pd{x)/3). □ 



< sup 

\ ze(B(x.2on(x)/3)) 



\uiz)\ Py{t > T,X, G {Bix,2pDix)/3))} + 
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4 Boundary Harnack inequality of A^'*^ on C^'^ ^ (C*^'^) open sets 

Next we assume that € dG and choose the coordinate system CSq. For x G M", r > 0, let 
A(x, o, r) be the box defined by 

A(x, a, r) = {y = (y, y„) S G : < - ^o{y) <a,\y -x\< r}. (4.1) 

The following result is a special case of Theorem 1.1. 

Lemma 4.1. Let 1 < a < /? < 2 and Zet G he a G^'^~'^ open set with characteristics rg = 1 
and A. Assume that k satisfies the conditions in Proposition \2.4\ Then there exist constants 
As = Agi^n, a, /3, A, Co, Ci, C2) < 1/2 and Ag = AQ{n, a, /?, A, Co, Ci, C2) such that 

<P,{X.^(„ ^^ G G} < A9/,(x)"-^ (4.2) 

for X G A(0,A8,^8) w'^^/i X = 0. 

Proof We assume that k = 1 because the proof is the same for the general case. Let p = 
(a - 1 + ((a + /? - 2) A l))/2 and define 

viiy) =Ua-i{y) +Up{y), 

where Ua-i and Up are functions defined in Proposition 12.41 Since p — a > (3 — 2, hy Proposition 
[231 there exists ki = ki{n,a, P,A) such that A(0, 2/ci, 2A;i) C B{0,ro) and 

A'^^\,{y)>0, y G A(0,A:i,fci). (4.3) 

Let (/> be a C^ function on C such that 

Hy) = \y\=yl + --- + yl-i, Iy| < i; i < </'(y) < 2, ly| > 1. 

Define 

V2{y) =ua~i{y) - Up{y)/{2Al) + Uk^^Al(t>{y). 

By Lemma 3.4 [12], we have {Aq!"^ (j){y)\ < A;2(p(y)^ " V 1), y G A(0, A;i, A;i), for some constant 
^2 = k2{n, a). Thus byp — a<l — a and Proposition 12.41 there exists m = m{n, a, /3, A) < ki/2 
such that 

A'^^\2iy)<0, yG A(0,m,/ci). (4.4) 

Since V2 > SAf on G \ A(0, oo, ki/2) and V2{y) < A4p{y)°'-^ for y G G n 5(0, ro) with y = 0, we 
have by applying (j3.7p and ()4.1ip 

^x{^r^(o,™,.,/2) e G \ A(0, oo, ki/2)} < 3-'Afpixr-' (4.5) 

for X G A(0, m, fci) with x = 0. 

Noticing that supyg^ l^i(y)l < and f i(y) > A^^p{y)"'^^ for y G G n S(0, ro) with y = 0, 
by ([3771) and ([Oil , we have 

^-{^-A(o.,n,.i) eG\A(0,m,A;i/2)} > 2-iAjV(a:)""\ x G A(0, m, fci) and £ = 0. (4.6) 

Combing (|4.5p and ([4.6p . we have 

^-{^rA(o,™,.,/2) e A(0,oo,^i/2) \ A(0,m,fci/2)} > e-^yl^ V(:e)°-i (4.7) 
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for X E A(0, m, ki) with x = 0. By (j4.7|) and (j3.12p . we can find a constant k2 = k2{ki,m, A) 
PAXr^,,^,,^/,, e A(0,fci,fci/2) \ A(0,m,fci/2)} > A;2^jV(^)"-i. (4.8) 
By (gS]), for X G A(0, m, /ci) with x = 

>P,{X,^(„_^ G A(0, 2A;i, fci) \ A(0, k,, k,)} 
>PAXr^,o,^,,^/,, G A(0, k,, k,/2) \ A(0, m, k,/2)}. 

sup PJX,^ G A(0,2A:i,fci)\A(0,A:i,A:i)} 

S/GA(0,fci,fci/2)\A(0,m,A:i/2) 
>k2k3Afp{xr-\ (4.9) 

where we use the fact that for some ks = k^{ki,m^ A) 

Py{Xr^,o,,,M) ^ A(0,2A:i,A;i) \A{0,ki,k^) } > k^, 

yeA{0,kuki/2) \A{0,m,ki/2). (4.10) 

which can be proved by the same calculation as (j3.25p . Setting Ag = ki, (j4.5p . (|4.9p and (|4.10p 
yield the first inequality of (j4.2p for x G A(0, fei, fci) with x = 0. 

Set ^3(2;) = V2{x)Lj.£G,\x\<i/2 + Ix€G,\x\>i/2i by Proposition l2.4l we can choose ki small enough 
such that 

A°/'7;3(y) < 0, yG A(0,fci,fci). (4.11) 
This estimate and Proposition 13.11 gives the second inequality of (j4.2p . □ 

Lemma 4.2. [Carleson estimate) Let 1 < a < /3 < 2 and /ei G he a C^'^~^ open set with 
characteristics tq = 1 and A. Assume that k satisfies the conditions in Proposition \2.4\ Let 
Q = G dG and assume that u > is a function on G which is not identical to zero, harmonic 
on GriB{Q, 1) and vanishes on dGriB{Q, 1) continuously. Then there exists a constant Aiq = 
Aio(n, a, /3, A, Co, Ci, C2) such that 

u{x) < Aiqu[xq), X £GnB{Q,l/2), (4.12) 

where xq = (0, 1/2) in the coordinate system GSq. 

Proof By chain arguments, we only need to prove (j4.12p for x G GriB{Q, 1/8). By multiplying 
a constant we can also assume that u{xo) = 1. Choose < 7 < a/{n + a) and define 

Bo = GnB{x,2p{x)), Bi = B{x,p{xy). 

Set 

B2 = B{xo, p(xo)/3), ^3 = B{xo, 2p(xo)/3) 

and 

n = inf{t >0:Xt^ Bo}, T2 = mf{t > : Xt ^ B2}. 
By (j4.2p and scaling, we can find a constant 5 = d{n, a, 13, A, Co, Ci, C2) such that 

Px{Xr^ edG)>6, X G C n B{Q, 1/4). (4.13) 

By Harnack inequality (j3.14p . there exists P' = (3{n, a, 13, Co, Ci, C2) such that 

u{x) < p{x)--^'u{xq), X G CnB(Q,l/4). (4.14) 
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Since u is harmonic on G fl B{Q, 1), we have for x E G n B{Q, 1/4) 

U{X) =E^{u{Xr,)IXr,eB,) + E^{u{Xr,)Ix^,iB,)- (4.15) 

We first prove that there exists constant /q > such that 

E^{u{Xr,)Ix^^iB,) < u{xo), xeG^n B{Q, 1/4). (4.16) 

Denote the Green function of (Xt) on an open set U by . For x G G'^/s l~l B{Q,l/4) 
satisfying 

\x - y\ < 2\z - y\, z G Bq, y ^ Bi, 
we have by Proposition 13.21 and the last conclusion in Lemma 13.41 

Ex{uiXri)Ix,_^^Br) 

A< \ I' I' K{z,y)G^o{x,z) 

=A[n, —a) / / j T— u[y)azay 

JBoJy€G,\y-x\>p(x)-y \z - y\^+'^ 

<2"+->^(n, -a) [ G^°(x,z) [ -S^I^^dzdy 

Jy<^G,\y-x\>pix)'< 



^JydG.. 



C2u{y) 

— uu 

a ^ 

u{y) 



\x - 



yeG,\y~x\>p(x)'i,\y-xo\>2p{xo)/Z \x - ?/|"+" 

flu 

X - 



dy 



u{y) 



J\y-xo\<2p{xo)/Z 

:=2"+2"c'2^^^(n, -a)p{xr{h + h)- (4.17) 



Similarly, 

u{xo) >E,^„{u{Xr^)Ix^^<f,B3) 

>2-(-+")Ci^(n,-a) / G^^(x,z) / ^^^^dzdy 

Jb2 JyeG,\y-x\>2pixo)/3 " ^1"+" 



>2-("+")Ci^6^(n,-a)(p(xo)/3)- / , ""^^l+j v- (4-18) 

JyeG,\y~x\>2p(xo)/3 ^0 - ^1"+" 

We have \y - x\ > 2-^p{x)^\y - xo\ if \y - x\ > p{x)'^ and x G 1/4). This and (fiJ8]) show 
that 

/i <2"+"/5(x)-^("+") ! ^^M}^dy 

Jy&G,\y-xo\>2p{xo)/?. FO " ^1"+" 

<22("+")3"^[i(Gi^(n, -a)p(xo)")-V(2;)"^^"+"^^i(2;o). (4.19) 
On the other hand, if p{x) < p{xq)/6, we have by Harnack inequality (j3.14p 



y-xo\<2p{xo)/3 F - y 



_u{xo) 

y-x\>p{xQ)/e 

<(27r)"2"+"yl7(p(^o)/6)-"n(xo). (4.20) 
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Combing (|4.17p - (|4.2U|) . we have for some constant c = c(n, a, (3, A, Cq, Ci, C2) 

E,{u{Xr,)Ix^^^B^) <cp(x)-^("+")7x(xo), X G G,(,„)/6 n B{Q, 1/4). (4.21) 

Noticing that a — 7(11 + a) > 0, by choosing Iq = lo{n, a, (3, A, Cq, Ci, C2) small enough, we get 
([TO]) from (|OTD . 

Suppose that there exists xi G GnB{Q, 1/8) such that ti(xi) > M = M(n, a, f3, A, Co, Ci, C2) 
> /q^' V(l + (5-^) (M will be fixed later). By KW\ . M > l^^' and u{xq) = 1 we have p{xi) < Iq. 
By (l4l^ . (f4T6]l and M > 1 + 

ii;,,KX,J/x.^eB,)>^M. 

From this inequality and (I4.13P we can find X2 G C such that 

|xi -X2I <p(xi)^, n(x2) > (1 - <52)-iM. 
Inductively, if S G fl B{Q, 1/4) for some k >2, we can find x^+i G C such that 

\xk+i - Xk\ < p{xk)\ u{xk+i) > (1 - 5'')-^u{xk) > (1 - 5'')-^M. (4.22) 

By (imil and (fO^j) . we have /)(xfc) < (1 - (52)''/'^'M-i/^'. Therefore, if (f02]) holds, we have 

fc-i 

\xk\ <\xi\ + Y, l^m - ^^1 < 1/8 + (1 - (1 - 5^f'P')-Hr^i^'. 

i=l 

Thus for M = (/q ^' V (1 + 6-^)) V (8^'(1 - (1 - J^)!//?')-/?')^ we can find x^ G C n B{Q, 1/4) 
satisfying K22h for all A; > 1. This gives a contradiction by noticing that lim^ — ^^y^ 7/(xfc) — 00 
and that u vanishes on dG fl B{Q, 1) continuously. Therefore supj^g^jp^j-g ^/s) ^(y) — 1^ 

Remark 4.1. Lei C and k satisfy the conditions in Theorem \l.ll then we can prove the hitting 
probability estimates in Lemma \4-l\ and the Carleson estimate in Lemma \4-S\ still hold. This 

— ana k = 1 



is due to that we have the same (super, sub) harmonic functions for k 
{see (jg-ip l and the term C'|x — y| does not destroy the (super, sub) harmonic functions which we 
construct above (c.f. Lemma \2.3\) . We omit the proof of this extension because it can be done 
by following the arguments for k & C^{G x G). Notice that function k in Theorem \l.l\ satisfies 
condition [Ci, C2, C3, — 1] for some constant Ci,C2,G^, and hence the Harnack inequality in 
Theorem \3.5\ holds. 



Before proving Theorem 1.1, we give some remarks on the assumptions of k and C. When 
■02 7^ 0, the condition of n in (jl.Sp is to study the reflected subordinate Brownian motion. 
However, due to the definition of the reflection point, we need C^'^ condition on C in Theorem 
1.1 when Tp2 7^ 0. Let C be a C^ open set in R". By the Appendix in [22], there exists 60 > 
such that, for any x G G^^^, there is an unique point ^(x) G dG satisfying |x — .^(x)| = p{x), 
C G C^(G^ ) and p eG'^{G'^ ). For x G G'^^, define the reflection point of x by 



x = 2^(x)-x. (4.23) 

When G is a C^'^ open set, and p are Lipschitz and C^'^ in a neighborhood of dG, respectively. 
The proof for the uniqueness of ^(x) is similar to [22]. The Lipschitz and C^'^ properties follow 
by the C^ case and the standard smooth approximation. 
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Proof of Theorem 1.1: First we assume that G is a C^'^~^ open set with characteristics 
ro < 1, A and k satisfies the conditions in Proposition 12.41 Let u > be a function on G which 
is not identical to zero, harmonic on G D B{Q,r) and vanishes continuously on dGnB{Q,r) for 
some < r < rg. By scaling and translation we can assume that r = 1 and Q = 0. Take the 
coordinate system GSq and denote 

Bo = AiO,As,As), Bi = A{0,2As,2As), r = tbo- 
By scaling, we can also assume that 

BiCB{0,l/3). (4.24) 

Write 

xo = (0,1/2), rci = (0,3yl8/2). 
By Harnack inequality (j3.14p . we have 

K^u{xo) <u{x) <kiu{xo), x£ A{0,2As,As)\A{0,As,As) (4.25) 

for some constant ki = /ci(n, a. A, Ci, C2, M). Next we assume that x G A{0, A^, Ag) with 
X = 0. Since u is harmonic on G fl B{Q, 1), we have by (|4.2|) and (|4.25|) 

u{x) =EMXr) > kY\{xo)P^{u{Xr) G A{0,2As,As)) 

>A^^k^\{xo)p{xr'\ (4.26) 

By the same calculation as (j4.18p . we have 

n(xi) >k2 [ I ""^^J y (4.27) 
Jy^B, \xi - yr+° 

for some k2 = A;2(n, q, Ci, C2, M). By (j4.24p and Proposition 13.21 we can also find a constant 
ks = k3{n,a,A) such that (c.f. (|3.18p ) 

E,{T)<k3P^{XreG\B^). (4.28) 

By definition of xi, we can find a constant ^4 = A;4(A) such that \z — y\ > k4,\xi — y\ for 2; G -Bo 
and y ^ Bi. Thus, by Proposition [321 Lemmas [331 Sll H^] and applying (^iM), (fOSl) . (|i:271l 
and (I4.28p . we have 

f Af \ C C K{z,y)G''^{x,z) 
u{x) =A{n, -a) / / . — — — u{y)dzdy 

J Bo JyeGnB-L \z — y\ 



^A( \ f f i^{z,y)G''°{x,z) 

+Ain,-a) / / -— — u{y)dzdy 

J Bo Jy^^GnBi \z - yr+" 



<C2^(n, -a)Aiou{xo)P^{Xr G B{) + G2A{n, -a)A:-("+") j [ ^-^^l^u{y)dzdy 

J Bo JyiBi \X1 - yr+" 

<C2^(n, -a)As,Aiou{x^)p{xr-^ + G^Ain, -a)k-^''+''^ E,{t) [ ^M^dy 

Jy^Bi kl - 2/1"+° 

<G2A{n, -a)u{xo)AQ{Aio + kik^^k3k-^''+''^)p{x)'''\ (4.29) 

Combing (j4.26p and (j4.29p we prove ()1.6p for x G A{0, Ag, Ag) with x = 0. By considering the 
coordinate system GSy for y G B{Q,2/3) fl dG, applying the arguments above and the Harnack 
inequality (j3.14p . we can prove (jl.6p . The general case can be proved similarly with the help of 
Remark 14.11 

□ 
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5 Boundary Harnack inequality of on Lipschitz domain 

To simplify notations, we assume that k = 1 in the arguments below because the estimates are 
the same for the general cases. Let G be a Lipschitz domain with characteristic tq and A, i.e., 
for each xq G dG, we can find a Lipschitz function T^q ■ M""^ — > M with Lipschitz coefficient 
not greater than A and an orthonormal coordinate system CS^q with which it holds that 

Gr\B{xo,ro) = {y = {yi,--- ,yn) : Vn > '["xoivi, ■ ■ ■ Vn-i) } B{xo,ro) . (5.1) 

The following hitting probability estimate is obvious for the Brownian motion case. Here we 
use capacity to give this estimate. We refer to [21] for more details about capacity and energy 
measure class 5*00 of symmetric Markov processes. 

Lemma 5.1. Let G be a Lipschitz domain with characteristic tq = 1 and A. Let xq = £ dG 

and choose the coordinate system GSq. Assume that A is a constant such that A(0,A,A) C 
Gn-B(0, 1). Then there exists a constant An = Aii{n,a,A) such that 

P^iXredG)>l/An, x = {0,xn), 0<Xn<A/2, (5.2) 
where r = inf{t > : Xt e A(0, A, A)"}. 

Proof By scaling, we may assume that A > 1/(3(1 + A)) without loss of generality. Denote the 
heat kernel of (Xt) by p{t,x,y). By Theorem 1.1 in |17j, there exists constant ki = ki{n,a,A) 
such that 

fei ii-''^" A j—^) < Pit, X, y) < fcr^ (t-"/" A _ [ ) , 0<t<L (5.3) 



Set F = A(0, A, AY n G and denote by {Yt) the killed process of [Xt) when hitting F. We know 
that the heat kernel p^{t,x,y) of (1^) is given by 

p'^{t,x,y) = p{t,x,y) - p{t - s,z,y)Px{X^ e dz,a e ds), (5.4) 

Jo Jf 

where a = inf{t > : Xt G F}. Noticing that A > 1/(3(1 + A)) and choosing 6 small enough 
we get by ([O]) and (fS^i]) 

/(t,x,y) > A:2(t""/"A < t < 1, x,y eG n B{0,6) (5.5) 



for some /c2 = k2{n,a, A). Set T = dG n B{0,6). Define the 1-potential kernel of (Yt) by 
Ui{x, y) = Jq°° e~^p^{t, X, y) dt and define for measure on G\F 

?7iV(x)= / U^{x,y)i^{dy). (5.6) 
Jg\f 

By Theorem 4.2.5 in [2T] and the continuity argument, there exists a 1-equilibrium measure 
supported on T such that 

U^M^) = Elie-''^), xeG\F, (5.7) 

where or = inf{t > : G F}. By problem 2.2.2 in [21], 

MT) = sup{/i(r) : G Soo, suppifi] C K, [7°^ < 1}- (5.8) 
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Direct calculations shows that /x = 6Irm{dx) G Sqq for 6 > 0. Choosing 6 small enough and 
applying the second inequality in (j5.3p . we get z^r(r) > k^in, a, A). Therefore by (jS.Sp and ([57 
we have for x G G fl B{0, 6) 



r JO 



>e "'^t'r(r) inf / po{t,x,y)dt > k4{n,a,A), (5-9) 

which implies that 

P^{ar < oc) > E^^ie-"^) > k4. (5.10) 

Noticing that Px{Xr G ^G) > P°(crr < oo), we get ([SJ]) for x G Gn B(0,(5). Thus we complete 
the proof by the Harnack inequality in Theorem 13. 5i □ 

Let ^ be a constant such that A(0, 2A, 2A) C G H B{0, 1) under the coordinate system CSq 
for a Lipschitz domain G with characteristic tq = 1 and A. Set 

Ko=A{0,A,A), Ki = A{0,2A,A)\Ko, K2 = G \ AiO,2A, A); (5.11) 
Hi ={Xr^, e ^1}, H2 = {X,,., G K2}. (5.12) 

Lemma 5.2. With notations defined in l[5.11\} and \5.1S^) . for any /c > 0, there exists a constant 
A12 = Ai2{n, a, A, A, k) such that 

Py{Hi) > Ai2p{yr\ lnp{y)f, y G Kq, \y\ < A/2. (5.13) 

Proof Let y G Kq with \y\ < A/2 and r = TB(y,p(y)/2)- We assume also that B{y, p{y) /2) C Kq. 
Otherwise (j5.13p can be verified by showing that Py{Hi) > ki{n,a,A,A). As the calculations 
in (|3.18p . we have by Lemma [ 



Py{Hi) >k2{n,a)Ey{T) [ [ dz 

>kk3{n,a,A,A)p{yr, (5.14) 

which gives ()5.13p for k = 0. Suppose (j5.13p holds for some k > 0. By the strong Markov 
property 



PyiH,)>Py{T<TK„ eK,: 

= [ P,{Hi)Py{Xr G dz) 
JKn 



Ko 

>h{n,a) [ P.{H,) ^'^]}^J z 

JzeKo\B(y,p(y)m \z-y\^+'^ 

. . ,N /■ p(y)°/o(2)"|lnp(2)|'= , 

>k5{n,a,A,A,k) / 1 mJa ^z. (5.15) 

JzeKo\B{y,p{y)m,\z\<A/2 \z - vl 

Direct calculation shows that Py{Hi) > A[2P{y)'^\^^ p{y)\^^^- Hence the proof completes by 
induction. □ 

Lemma 5.3. With notations defined in k5.11]) and Ii5.12\) . there exists a constant A13 = 
Ai3(n, a, A) such that 

PyiH2) < AisPyiH,), y£Ko, y = 0. (5.16) 
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Proof To simplify the arguments, we assume that G is a special Lipschitz domain. By scaling, 
we assume that ^ = 1 in (jS.lip . For i > 1, set 

J, = A(0, 2-\n) \ A(0, 2-^-\n), = i - ^ ^ 1 

and ro = ri. Define for i >1 

di = sup P,{H2)/P,{H^), Ji = A(0,2-%r,_i), Ti = rj.. (5.17) 

z€.Ji 

By Harnack inequality, each di is finite. Noticing that Tj < t^q and applying the strong Markov 
property, we have for z G J, and i > 2 

i— 1 « 

< / P,{Xr, G dw)P^{H2) + P,{Xr,eG\ UiL\jk) 

k=l ■'■^k 
i—1 „ 

<Y.^k Pz{Xr, G dw)PUHl) + Pz{ Xr,£G\ UlJ^Jk) 
k=l -^-^k 

<( sup dk)Pz{Hi) + P,{Xr, G G \ U^-\ Jfc). (5.18) 

i<fe<i-i 

Define uq = 0, cti = inf{t > : \Xt — Xq\ > 2^*} and define by induction cJm+i = cti o ^^.^ for 
m > 1. Similar to the calculation of (j3.18p . we have for some constant ki < 1 independent of i 
such that ^ ^ 

Pyj{X„^ G Ji) < 1 - Pw{Xa^ G U*^~\ Jfc) < tt; G Jj. 

Therefore, for z € Ji and positive integer /, we have by the strong Markov property for 
Pz{r^ > au) <Pz{X^, €Ji,l<k<li) 

= f Pz{Xa, eJi,l<k<li-2, X,^^_^ G dw)Py,{X„, G Ji) 
<Pz{X„^ e Ji,l <k <U - l)ki < k[\ (5.19) 

Recall that x is the first n—1 coordinate of x. On {Xr^ G G\U^"L^^ J^, Tj < an} with Xq = z E Ji, 
we have l^o-^ — ^o-ol > 5^ ~ 2~* for some 1 < k < li which implies for some 1 < k' < li 

Therefore, we have for some io ^ 2 

{Xr^ eG\u{z}^Jk, Ti < au} 
CuLi -X,^_J > 1/(100/^3), G J^}, i > io, (5.20) 



and hence 



P,{Xr^ e G\U^Jfc, Ti < an 

li 



<^P,(|X., -X.,_J > 1/{100H'),X,^_^ G Ji) 



k=l 
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<li sup P,{\X^^ \ > 1/(100/^3)) 

<k2li2-'^'{li^T. (5.21) 

The proof of the last inequahty above is similar to (j4.17p . By (j5.19p . (j5.2ip and choosing / big 
enough, we have for z £ Ji and i > io, 

P,{Xr, G G \ U^-\ Jfc) < + k2m~''\li^r < A:32-"'«=^"+^ (5.22) 

By (f5J8D . (102]) and Lemma[01 for z e Ji and i > io 

Pz{H2)/P,{Hi) < sup 4 + Pz{Xr, G G \ yff72iJk)/Pz{Hi) < sup 4 + ki/i''- 

l<fc<j-l l<A:<j-l 

This implies that 

i 

di < sup 4 + ^4 2^-*-/^^— 4 + 3A;4, 

l<fc<io-l ^ l<fc<fio-l 

which completes the proof of this lemma. □ 

Remark 5.1. One may use the method in 110^ to give a better estimate of i5.13\) . The proof of 
Lemma 15.31 is an adaption of the Brownian motion case. 

Proof of Theorem 1.2: In the proof of Lemma 14.2^ we only use the G^'^~^ property in 
(j4.13p . Thus we can prove the Carleson estimate for the Lipschitz case with Lemma [5. II in place 
of M.lSp . Therefore, we can prove Theorem 1.2 by the standard arguments of BHI with the help 
of Theorem 13.51 and Lemma 15.31 □ 

6 Boundary Harnack inequality of A"/^ 

When G = W, A^''^ is the fractional Laplacian A"/^. Recall that Wp{y) = for y G M!^. We 

extend these functions to M" by taking zero on M" \ . Next we give a formula of A'^/^tUp. 
Integration by parts formula shows that for < p < a 

yV — rP r°° ip - 1 

\y-x\^+^ Jo |y-l|i+" 

= limxP-" ( [ ' , dy + r - y^-ll-dy] 

Thus, for n = 1 we have 

r\ <y. — p — 1 'p- — 1 

^''I'^Wpix) = A{1, -a)^j;P-" / — ^| dy, < p < a. (6.2) 

Jo |y — ir 

Applying spherical coordinate transform from (yi, • • • , y„) G M*^ to (r,9i, - ■ ■ , On-i) G [0, oo) x 
[0,7r]"~^ X [0,27r), this gives for n > 1 

A{n,-a)~^A'^/'^Wp{x) 

= lim/ ^^^#^^^2/ 

"^iO J\y~x\>e P y\ 
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: lim d^i • • • de„.-2 [ " ^{Oi, ■ ■ ■ en-2) de. 

Jo Jo Jo 

I H>o Icosei ) 



j/lt--^2!_i>p\ cos^ t^i— ; — dt 

ll'^ cosSl^^J f \l + a 

-oo I'' COS 6*1 I 

- / ^ / dOi--- den-2 / ^(^1, • • • en-2) (cOs" 

a Jo |y-ir Jo Jo Jo 



"Jo |y-lr J\y\=l,yn>0 
where (p{6i, ■ ■ ■ On— 2 

) = sin'^ ^ 6*1 sin" ^ 6*2 • • • sin 

6n—2j iTi{dy) is the (n — l)-dimensional Lebesgue 
measure and we use the following transform in the calculation above 

r = t — x„/ cos ^1, 0iG[O,7r/2); — r = t + x„/ cos 6*1, ^1 G (7r/2, vr]. 

Denote for < p < a and n > 1 

A(n,a,p = / . — dy / y„ m{dy), 

a Jo |y-lr J\y\=i,yn>0 



Tt \ A{n,-a) ( f 

A{n,a,p) il+p 

a V Jo 



-dy ] I y^ m{dy) 



y|=i,s/n>o 



with convention that m{dy) is the Dirac measure for n = 1. By (j6.3p . we have the following 
Lemma. We notice that the case p = a/2 below has been obtained in Example 3.2 of [7]. 

Lemma 6.1. Let < p < a < 2, we have 

A"/2u;p =A(n, a,p)wp-a, x G Ml, (6.4) 

A'^,^Wp =A{n,a,p)wp-a, x G M" . (6.5) 

Formula (j6.5p is another version of (j3.ip for < p < a. By Lemma [6. II and (j3.ip we see that 

A"'/'^Wp<0, -l<p<a/2; A^/^y^p = 0, p = a/2; A^/^^^j^ > 0, a/2<p<a. (6.6) 

Let K be a symmetric function on M" x M" such that 

Ri < K{x,y) < R2, \K{x,y) — k{x, x)\ < R-^lx — y\, x,y gM."^ (6.7) 

for some constants Ri,R2,R^ > 0. In what follows, notation (Xt) is for the stable-like process 
on M" associated with A 2'''^. Harmonic functions of (Xt) is again defined by (|2.ip . 

Lemma 6.2. Let < a < iVa < /3 < 2 and D a C^'^~^ open set in M" with characteristics rg = 
1 and A. Let k be a symmetric function on M" x M" satisfying 1^6. 7| ). T/ien /or a/2 < p < a and 
Q G i/iere exist function Up and constants A13 = Ai^{A), An = Ai4{n, a, /3, A,p, Ri, i?2) -R3) 
SMc/i that 

AiilDnB{Q,2/3)p{xy < Up{^) ^ '^nlDnB{Q,2/3)p{xT, X G M", (6.8) 
At'«-up(x) > Ai4p(x)P-", X G i:>nS(g,l/^i4), a/2<p<a, (6.9) 



Ai4p(x)^-"/2-i, X G Z)nB(g,l/2), /3<l + a/2, 
|At''^u„/2(x)| < <( ^i4|log/9(x)|, X G Z)nB(Q,l/2), /3 = l + a/2, (6.10) 
^14, X G £>nB(Q,l/2), /3>l + a/2. 
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Proof Following the calculations in Lemma 12.11 Lemma 12.21 and Proposition 12.41 we can prove 
this lemma with the help of ()6.4p and ()6.6p . We omit the details of the proof because, by noticing 
that Wp = on M" \ , the calculation is essentially on D which is the same as the regional 
fractional Laplacian case. □ 

By Lemma 16.21 and following the arguments in Proposition 14.11 and Theorem 1.1, we can 
prove the following results. 

Lemma 6.3. Let 0<a<l\/a<(3<2 and D a C^'^~^ open set in M" with char- 
acteristics ro = 1 and A. Assume that Q = £ dD and k is a symmetric function on 
M" X satisfying \6. 7p . Then there exist constants Ai^ = ^15(71, a, /3, A, i?2, -Rs) < 1/2 
and AiQ = ^i6(n, a, /3, A, Ri,R2, R3) such that 

A^ipixr/' <P.{^.,,o,.,5.^i5) ^ A(0,2Ai5,^i5)} 

<^-{^-A,0,.,5.^15) ^ ^} ^ ^16P(X)"/^ (6.11) 

for X S A{Q, Ai5, A15) with x = under CSq. 

Theorem 6.4. Assume that a, (3, D and k satisfy the same conditions as in Lemma \ 6.3l 
Let Q E dD and r € (0, tq). Assume that u > is a function on D which is not identical 
to zero, harmonic on D D B{Q,r) and vanishes on D'^ n B{Q,r). Then there exists constant 
C = C{n, a, A, Ri, R2, R3) such that 

44 ^^4^' forx,yeDnB{Q,r/2). (6.12) 

Remark 6.1. By taking G = D, all the conclusions in Section 3 can be extended to A^''^ in a 
similar way, where the reflected stable process is replaced by the stable-like process. The Carleson 
estimate for A^''^ can be proved by the same method as in Lemma \4-^ We remark that to prove 
the boundary Harnack principle of A 2''^ on open sets, we need the method in [8] to get the 
Carleson estimate, where the explicit Poisson kernel can be replaced by the sharp estimates as in 
Theorem \6.4\ can be generalized to operator Aq when we further assume that D C D G G. 
The proof of this generalization is the same as the case G = M"" except that the constant depends 
also on the distance between D and dG. 

Remark 6.2. Since {{wp)p^i, (wp)p>i) wi is the (super, sub )harmonic function of Laplacian 
on half spaces, by the Harnack inequality in ]30^ and the method in this paper, we can prove the 
explicit BHIfor A + A"/^ on C^'^ open sets which gives p{x) order decay for harmonic functions 
near the boundary. With the help of this fact we can prove that the Green functions 0/ A + A"/^ 
and A are comparable on a G^'^ bounded open set. 
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